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Abstract 

. We present the results of a study of the semileptonic decays of the Bc-meson and 

I the lowest lying doubly heavy baryons using the relativistic quark model. We do 

■ not employ a heavy quark mass expansion but keep the masses of the heavy quarks 

. and hadronss finite. We calculate all relevant form factors and decay rates. 

: 

'"J ■ The semileptonic decays of heavy mesons and baryons are ideally suited to extract the 
, Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The heaviest flavored bottom- 
1 charm i?c-meson was observed by the CDF Collaboration [||] in the analysis of the decay 
mode Be — >■ J/tplv. The discovery of the i?c-meson raises hopes that doubly heavy 
flavored baryons will also be discovered in the near future. The theoretical treatment of 
the systems with two heavy quarks is complicated by the fact that one cannot make use 
^ , of an expansion in terms of the inverse heavy quark masses. Previously, nonrelativistic 
O ' potential models, diquark approximation, QCD sum rules and nonrelativistic QCD have 
^ \ been used to describe the spectroscopy of doubly heavy baryons and to estimate the 
' inclusive and some exclusive decay modes of such systems (for review, see and 
p • references therein). 

We present here our recent results |^ of exploration of the semileptonic decays of the 
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D doubly heavy i?c-meson and the lowest lying doubly heavy S-baryons within a relativistic 
constituent quark model. This model P] can be viewed as an effective quantum field 
' theory approach based on an interaction Lagrangian of hadrons interacting with their 
^ constituent quarks. Universal and reliable predictions for exclusive processes involving 
both mesons composed from a quark and antiquark and baryons composed from three 
quarks result from this approach. The coupling strength of hadrons H to their constituent 
quarks is determined by the compositeness condition Zh = |^, O] where Zh is the wave 
function renormalization constant of the hadron. The quantity Z^ is the matrix element 
between a physical particle state and the corresponding bare state. The compositeness 
condition Zh = enables us to represent a bound state by introducing a hadronic field 
interacting with its constituents so that the renormalization factor is equal to zero. This 
does not mean that we can solve the QCD bound state equations but we are able to 
show that the condition Zh = provides an effective and self-consistent way to describe 
the coupling of the particle to its constituents. One starts with an effective interaction 
Lagrangian written down in terms of quark and hadron variables. Then, by using FejTiman 
rules, the S-matrix elements describing hadron-hadron interactions are given in terms of 
a set of quark diagrams. In particular, the compositeness condition enables one to avoid 
the double counting of quark and hadron degrees of freedom. This approach is self- 
consistent and all calculations of physical observables are straightforward. There is a 
small set of model parameters: the values of the constituent quark masses and the scale 
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parameters that define the size of the distribution of the constituent quarks inside a given 
hadron. The shapes of the vertex functions and the quark propagators can in principle 
be determined from an analysis of the Bethe-Salpeter (Fadde'ev) and Dyson-Schwinger 
equations, respectively, as done e.g. in In the present paper we, however, 

choose a more phenomenological approach where the vertex functions are modelled by 
Gaussian forms and the quark propagators are given by local representations. We have 
demonstrated in our papers [|T3] that the relativistic constituent model is consistent with 
the heavy quark symmetry in the limit of infinite quark masses. 

We start with the effective interaction Lagrangian which describes the coupling be- 
tween hadrons and their constituent quarks. For example, the couplings of the bottom- 
charm hadrons into their constituents are given by 



Ant(a;) = gB, Bc{x) Jb,{x) + [g~,^ ^bc{x) + h.c.) , (1) 

Js,S^) = jdxi j dx2j dx3^s,S^; xi, X2, X3) {Ca2{x2) Qa^ixs)) e'''''^''', 

where q = u,d. The vertex function is taken to be invariant under the translation 
X — > x + a which guarantees Lorentz invariance for the interaction Lagrangian Eq. (|l]). 
The matrix elements of the weak current bO'^c are written in our approach as 
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where K = kf + {ki + k2) + k^. The quark propagator is chosen to have a local form 

^iif^) = — —Ti {i = u,d,s,c,b) (4) 

TTij— /C 

with rrii being a constituent quark mass. We assume that uih < Y17=i "^g, (n=2 for mesons 
and n=3 for baryons) in order to avoid the appearance of imaginary parts in the physical 
amplitudes. This is a reliable approximation for the heavy pseudoscalar mesons. The 
above condition is not always met for heavy vector mesons. We shall therefore employ 
equal masses for the heavy pseudoscalar and vector mesons in our form factor calculations 
but use physical masses for the phase space. 

Generally, is a function of extremal momenta too, however, in the impulse ap- 
proximation employed in our approach, we assume that it only depends on the sum of 
relative momentum squared. 

The coupling constants gn are determined by the so called compositeness condition 



proposed in [^] and extensively used in |]T0| . The compositeness condition means that the 
renormalization constant of the meson field is equal to zero Zh = 1 — = 0, where U'fj 
is the derivative of the mass function. The meson-mass functions are defined as 



np(p2) 
ny(p^) 



{27r)H 



9' 



-F)tr 



(5) 
(6) 



In the baryon case the compositeness condition may be rewritten in a form suitable for 
the determination of the coupling constants: 



-12^, 
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7"7' (^1 + P) Y ^91 + ph^^' tr {S,, (k, + k2)r ^93 (^2) 1'- 



(9192^3) = {bcq) , {bsq) , (esq) 



D'icq = 7°7' ^9(^2) 7^7' tr (7" S,{k^ + p) S,{k, + p) 7^ 5e(A;i + ^2)) . 

The calculational techniques are outlined in Refs. |^, |^. 

Before presenting our numerical results we need to specify our values for the constituent 
quark masses and shapes of the vertex functions. As concerns the vertex functions, we 
found a good description of various physical quantities @, |^, adopting a Gaussian form. 
Here we apply the same procedure using ^H{k\) = exp{— A;|;/A|^} in the Euclidean region. 
The magnitude of characterizes the size of the vertex function and is an adjustable 
parameter in our model. The Ah parameters in the meson sector were determined by 
a least-squares fit to experimental data and lattice determinations. The quality of the fit 
may be seen from Table 1 for the leptonic decay constants. The nucleon A^v parameter was 
determined from the best description of the electromagnetic properties of the nucleon . 
The Afj parameters for baryons with one heavy quark (bottom or charm) are determined 
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by analyzing available experimental data on bottom and charm baryon decays . Since 



there is no experimental information on the properties of doubly heavy baryons yet we 
use the simple observation that the magnitude of is increasing with the mass value of 
the hadron whose shape it determines. Keeping in mind that Ajy = 1.25 GeV, Aq^^ = 1.8 
GeV and A^^ = 2.43 GeV, we simply choose the value of Aqq^ = 2.5 GeV for the time 
being. We found that variations of this value by 10 % does not much affect the values of 
form factors. We employ the same values for the quark masses (see, Eq.(^) as have been 
used previously for the description of light and heavy mesons (baryons) |P, |l^. Note that 
the values for the light quark masses in the meson case are different than in the baryon 
case as a result of the lack of confinement in our approach. 
We thus use (in GeV) 

ruu rris rric rrib 



(0.235) 0.420 (0.333) 0.570 1.67 5.06 



Att A/< Ap Ap^ Aj/^ Ab Ab, Ab, At . . 
1.16 1.82 1.87 1.95 2.12 2.16 2.27 2.43 4.425 ^ ' 

The resulting form factors are approximated by the interpolating form 

It is interesting that for most of the form factors the numerical fit values of Oi and 02 
obtained from the interpolating form ([TD|) are such that the form factors can be represented 
by monopole function in the case of i?c-nieson and dipole formula in the case of doubly 
heavy baryons: 

f{<f) ^ ^f] , (n = 0,l). (11) 

The values of my in this representation are very close to the values of the appropriate 
lower-lying [qq') vector mesons (m/)*=2.11 GeV for (c-s)-transitions and rriB* ^ mB^=QA 
GeV for (b-c)-transitions). In Fig.l we show two representative form factors and their 
dipole approximations. It is gratifying to see that our relativistic quark model with the 
Gaussian vertex function and free quark propagators reproduces the monopole in the 
meson case and the dipole in the baryon case for most of the form factors. 

Finally, in Table 2 we present our predictions for the branching ratios of the semilep- 
tonic i?c-decay rates and compare them with other approaches. In Table 3 the predictions 
for the decay widths of the doubly heavy S-baryons are given. We compare them with the 
free quark decay widths. One notes that the rates for the exclusive modes "^cc + 

and Sfoc '^bs + are rather small when compared to the total semileptonic inclusive 
rate estimated. The remaining part of the inclusive rate would have to be filled in by 
decays into excited or multi-body baryonic states. Note that the smallness of the exclu- 
sive/inclusive ratio of the above exclusive modes markedly differs from that of the mesonic 
semileptonic 6 — > c transitions, where the exclusive transitions to the ground state S-wave 
mesons B ^ D, D* make up approximately 66% of the total semileptonic B ^ rate 
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2T| . For Afe — i> Ac transitions one expects even higher semileptonic exclusive- inclusive 
ratio of amount 80% Note that the rate for — ^ + li^ is of the same order 



of magnitude as the rates calculated for the corresponding double heavy mesonic decays 
Be Tic + W and Be J/"^ + 11^ The QCD sum rule and potential model predictions 
for the rates of — > Sec + and Sf,c — > Scs + W given in exceed our rate predictions 
by factors of 10 and 3 respectively. In fact, the exclusive semileptonic rates given in 
tend to saturate the inclusive semileptonic rates. 

In Table 4 we present values for the invariant form factors at g^j^ = and g^^^^^ = 
(mj — mj)^. Note that the values of the axial vector form factor are rather small for 
the two decays S^c — >■ Sec + and S^c '^bs + ^J^- This provides for a partial explanation 
of why the rates of these two modes are small compared to the inclusive semileptonic 
rate. Also the zero recoil values of the vector form factors are significantly below the 
value of one which one would expect from a naive application of the heavy quark limit. 
The smallness of the vector form factors provide for the remaining explanation of the 
smallness of the predicted respective rates. We mention that the QCD sum rule and 
potential model estimates of the zero recoil values of both the vector and axial vector 
form factors and F-^ given in are close to one. In the model of the form factors 
F2 and are set to zero. In our approach we find that the numerical values of F2 and 
F^ are quite small compared to those of F^ and F^ in all cases when expressed in terms 
of the mass scale [rrii + mj). 
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Table 1: Leptonic decay constants fn (MeV) used in the least-squares fit. 



Meson 


This model 


Expt. /Lattice 


n+ 

K+ 

D+ 

Dt 
B+ 

Bt 
Be 


131 
160 
191 
206 
172 
196 
360 


130.7 ±0.1 ±0.36 

159.8 ± 1.4 ±0.44 

19lli 
206+^^ 
172111 

171 ± 10tlt%^' 
360 


J/^ 
T 


404 
711 


405± 17 
710±37 



Table 2: Branching ratios BR(%) for the semileptonic decays Bt Hl'^u, calculated 
with the CDF central value t{Bc) = 0.46 ps 0. 



H 


This model 


115 


|16| 


|ic| 


m 


17 


|18| 


Vceu 


0.98 


0.8±0.1 


0.78 


1.0 


0.15(0.5) 


0.6 


0.52 




0.27 














J /ip ev 


2.30 


2.1±0.4 


2.11 


2.4 


1.5(3.3) 


1.2 


1.47 


J ji\}rv 


0.59 














D^eu 


0.018 




0.003 


0.006 


0.0003(0.002) 






D^TU 


0.0094 
















0.034 
0.019 




0.013 


0.019 


0.008(0.03) 






B^eu 


0.15 




0.08 


0.16 


0.06(0.07) 






B*^eu 


0.16 




0.25 


0.23 


0.19(0.22) 






B'^eu 


2.00 


4.0 


1.0 


1.86 


0.8(0.9) 


1.0 


0.94 


Bf e V 


2.6 


5.0 


3.52 


3.07 


2.3(2.5) 




1.44 
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Table 3: Calculated decay widths of lowest lying = 1/2+ doubly heavy S-baryons. 
Inclusive widths are calculated using the current quark pole masses. 



Mode 


Decay widths, ps ^ 


RTQM 


Inclusive width 


" be ^ '— ' cc ~l~ ^ ^ 


0.012 


2-ro(6^c)=0.162 


be — ^ bs + ^ ^ 


0.043 


Tq{c^ s)= 0.122 


' cc ^ '— ' cs ~t~ ^ ^ 


0.224 


2-ro(c^s)=0.244 



Table 4: Values of FY and F^ form factors at maximum and zero recoil. 





^bc ^ "CC 


^bc ^ ' — 'bs 


"CC "CS 




FY F^ 


Fi^ Fi^ 




0.46 -0.091 


0.39 0.061 


0.47 0.61 




0.83 -0.086 


0.58 0.065 


0.59 0.77 
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Figure 1: Upper panel: Form factor {q^) (solid-dotted line) for he — > hs transitions 
and its dipole approximation (solid line) with rrics — 2.88 GeV; Lower panel: Form factor 
{q^) (solid-dotted line) for he — > cc transitions and its dipole approximation (solid line) 
with mfec = 6.81 GeV. 
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